Magnetic-field symmetries of mesoscopic nonlinear conductance 
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We examine contributions to the dc-current of mesoscopic samples which are non-linear in applied 
voltage. In the presence of a magnetic field, the current can be decomposed into components which 
are odd (antisymmetric) and even (symmetric) under flux reversal. For a two-terminal chaotic 
cavity, these components turn out to be very sensitive to the strength of the Coulomb interaction 
and the asymmetry of the contact conductances. For both two- and multi-terminal quantum dots 
we discuss correlations of current non-linearity in voltage measured at different magnetic fields and 
temperatures. 

PACS numbers: 73.23.-b, 73.21. La, 73.50.Fq 
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I. INTRODUCTION 

Symmetries are of fundamental interest in all fields 
of physics. In linear irreversible transport the Onsager- 
Casimir 1 symmetry relations are important since they 
relate different transport coefficients. Here we are con- 
cerned with electron transport on the mesoscopic scale. 
This scale emerges when the distance carriers can travel 
without losing their phase coherence becomes compara- 
ble to the dimensions of the sample. It has long been 
understood that the Onsager-Casimir relations are not 
restricted to the macroscopic domain but extend to lin- 
ear transport on the mesoscopic scale. For the lin- 
ear conductance G a /3 = dI a /dV/3 Onsager-Casimir im- 
plies that under field reversal we have the symmetry 
G a p{&) = Gp a (—&). An important condition is that 
voltages are measured at contacts which are sufficiently 
large such that they can effectively be considered as equi- 
librium electron reservoirs^. 

Since symmetries are important it is crucial to under- 
stand their limit of validity. Are Onsager-Casimir rela- 
tions strictly valid only in the linear transport regime? 
What could cause their breakdown? Is it possible to 
quantify and measure the departure from symmetry? In 
this work we report recent theoretical and experimental 
progress on these questions using as an example electron 
transport through a chaotic quantum dot. We extend 
earlier discussions to correlations of current non-linearity 
in voltage measured at different magnetic fields and tem- 
peratures. Interestingly, as we will now discuss, these 
questions are related to the role Coulomb interactions 
play in non- linear transport. 

In the mesoscopic regime the linear conductance 
depends on quantum interference^ and depends on 
Coulomb interactions. A change in the Hartree potential 
is similar to a small change in the shape of the sample or 
in its impurity configuration and thus needs no separate 
discussion. In contrast Hartree-Fock terms^ can mod- 
ify the conductance significantly. A well known example 
of such an interaction effect is the physics of Coulomb 
blockade which becomes relevant at low temperatures if 
the contacts are pinched-off. However, if the conduc- 
tance of the quantum dot with ballistic contacts is large, 




FIG. 1: Quantum dot with magnetic flux $ and dc bias volt- 
ages Vi,2 at the contacts and Vo at the gates with capacitance 
C. 



G = (e 2 /h)g ^> e 2 /h, Fock terms give only a small rel- 
ative correction < 1/g to sample-specific quantum fluc- 
tuations. In such open dots weak localization correc- 
tions (WL) or universal conductance fluctuations (UCF) 
remain universal at low temperatures. In particular, 
Coulomb interactions do not affect WL and just slightly 
(< 1/g) modify UCF only at elevated temperatures^. 
Since these effects due to many-body physics in open 
quantum dots appear only as small corrections to non- 
interacting theory, one could conclude that Coulomb in- 
teraction effects are unimportant. 

However, Coulomb interactions play a much more im- 
portant role in non- linear transport. Experimentally 
transport is characterized by the full conductance I/V 
rather then the differential conductance, dl/dV ^ I/V. 
In particular, here we are interested in (l/2)d 2 I/dV 2 
(at V = 0) which we call the second-order non-linear 
conductance or simply the non-linear conductance. In 
this regime it is the Coulomb interaction that breaks the 
symmetry of the Onsager-Casismir relations. Conversely, 
measurement of transport non-linearity provides a tool to 
determine interactions. 

Initial discussions of non-linear dc- transport^ 7 - and 
the rectification of an ac-applied voltage^ in mesoscopic 
samples were addressed without taking interaction effects 
into account. Interestingly, inclusion of Coulomb interac- 
tion does not only modify these phenomena but also leads 
to a qualitatively new effect: namely, even in a two ter- 
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FIG. 2: Normalized fluctuation of the (anti)symmetric nonlin- 
ear conductance (Q a )Gs as a function of flux "J>/$ c ,$ c <C $o 
for strongly interacting dot with symmetric contacts, fromii. 



minal conductor it is possible to have a component to the 
current which is odd in magnetic field /($) 7^ /(—<&). In 
generic mesoscopic conductors, in chaotic cavities (quan- 
tum dots) and in metallic diffusive conductors this effect 
vanishes when an average is taken over many samples 
of slightly different geometry or impurity configuration 
(mesoscopic averaging)2ii£. Thus in generic mesoscopic 
conductors this is purely a quantum effect. 

Sanches and Biittiker investigated the effect of inter- 
actions on current-voltage characteristic of open chaotic 
(ballistic or diffusive) cavities, see Fig. 1, in the pres- 
ence of a magnetic flux $ comparable to a flux quantum 
$0 = eh/c£ They found the dependence of the fluctua- 
tion of the antisymmetric component on the interaction 
strength and the dependence on the number of ballistic 
channels transmitted through the contacts of the cav- 
ity. Spivak and Zyusin explored this asymmetry at small 
fields, $ <C $0, and weak interactions in open diffusive 
samples 10 . Under these conditions the asymmetric fluc- 
tuations are linear in flux and interaction strength. Po- 
lianski and Biittiker— present a theory which describes 
the entire crossover from the regime of low magnetic fields 
to the regime considered in Ref. H in which the flux 
through the sample is comparable to a flux quantum $0. 

The experiments measured field asymmetry in a wide 
range of structures: in nanotubes 12 , ballistic billiards 1 ^, 
ballistic quantum dotsii, ballistic^ and diffusive 1 ^ 
Aharonov-Bohm rings. This motivated further research 
on non-linear quantu m 17 i 18 and classical^ effects. The 
experiment by Zumbiihl et al*^ investigated this asym- 
metry in chaotic ballistic dots at very low voltages. The 
experiment demonstrated that the asymmetry vanishes 
on average and depends linearly on flux at sufficiently 
small fields. The experiment also investigates the depen- 
dence of the asymmetry on the number of channels in the 
contacts of the sample. Some experimental results, espe- 
cially, the magnitude of the asymmetry, if four or more 



channels are open, agree well with existing theories ^ 10 ' 11 . 
A more detailed investigation of the dependence on var- 
ious external parameters is still necessary. 

Theoretically Ref. [ll| considers fluctuations of both 
anti-symmetric and symmetric component of non-linear 
conductance 
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through a 2-terminal quantum dot with arbitrary inter- 
action strength, flux, temperature and dephasing— . The 
fluctuations (square root of the variance) of these two 
quantities are shown in Fig. [5] as a function of the flux 
applied to the cavity. Interestingly, similarly to linear 
transport (see references in 2 ^), the crossover of Q a and Q s 
from low to high magnetic fields is shown to depend on 
the flux <f> c <C $0 and not on the flux quantum $0 char- 
acteristic for open diffusive samples. Since the dot is es- 
sentially zero-dimensional, the dwell time an electron 
spends inside is much larger then the ergodic time r erg 
necessary to explore the phase space of the dot. There- 
fore an electronic trajectory gains a flux ~ $0 during 
Td/Teig S> 1 random explorations of the dot already at 
$ &o/ \Z T d/T C r g <C $o- Fluctuations of Q a , Q s occur on 
the background of an ensemble averaged anti-symettrized 
(Qa) = and symmetrized (Q s ) = non-linear conduc- 
tance. The magnitude of these fluctuations as a function 
of flux $/$ c is shown in Fig. [2] for the limit of strongly 
interacting dot with symmetric contacts (corresponding 
to the experiment^). 

However, if either the Coulomb interaction is not very 
strong or ballistic contacts are unequal, Ni ^ N 2 on Fig. 
HI the field-asymmetry of the non-linear signal may be 
strongly reduced compared to strongly interacting sym- 
metric set-up. In practice, to observe a strong violation of 
Onsager relation o 12 ' 13 ' 14 ' 15 ' 16 one might want to consider 
not the amplitude of Q a , but rather the relative asym- 
metry in the signal, denoted here A = S a /Gs- Below we 
discuss the optimal regime to maximize this ratio. We 
find how magnetic field and temperature affect fluctua- 
tions of this asymmetry parameter and find the contact 
asymmetry that maximizes A. For sufficiently strong in- 
teraction, symmetric contacts turn out to be the optimal 
regime. Surprisingly, however, if the interactions are not 
very strong an asymmetric set-up, N\ ^ N%, is more ad- 
vantageous. 

Since both field-components G a ,Gs were recently ex- 
plored experimentall y 14 ' 16 , we consider the statistics of 
Q a and Q s and their correlations for arbitrary interac- 
tions. Particularly, in the non-interacting limit we find 
an universal relation between UCF and fluctuations of Q s , 
independent of temperature and applied field. In view of 
a multi-terminal experiment^, we generalize our treat- 
ment of an asymmetric component Q a to a multi-terminal 
dot and consider non-linearity not only with respect to 
a contact, but also to the gate voltage. In the end we 
discuss interaction effects in linear vs. non-linear trans- 
port, applicability of existing theories to experimental 
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data and provide a formula for fluctuations of the asym- 
metry A in the experimental regime of Ref. Q 

The paper is organized as follows: Sec. |TT]describes the 
model and the basic steps of our approach. We present 
our results in Sec. IHII and discuss them in Sec. IIVI 



II. MODEL 

The 2D chaotic quantum dot, see Fig. [TJ is biased with 
dc voltages V a , a — 1, M at contacts with N a ballistic 
channels, and by the voltage V at the gates with capaci- 
tance C . Chaos in the dot is either due to diffusive scat- 
tering off impurities, when the mean free path is smaller 
then the size of the dot, I <C L, or due to the irregular 
shape of a ballistic dot with chaotic classical dynamics, 
I ^> L. The dot is in the universal regim e 20 ! 21 , when the 
Thouless energy Et = h/r CIg is the largest energy-scale 
in the problem. Typically £rh ~ fi min {Ivf/L 2 , vf/L}, 
see the review^ - for details. The mean level spacing (per 
spin direction) A = 2ith 2 /(rn* Area) and the number 
of conducting channels N together define the dwell time 
r d = h/(NA) ^> r crg which an electron typically spends 
inside the dot. 

Scattering is assumed spin-independent and this spin 
degeneracy is explicitly accounted for by the factor v s . 
The number N of ballistic (orbital) channels in the con- 
tacts is assumed to be large, N ^> 1, so that the prob- 
lem can be considered analytically. We use Random Ma- 
trix Theory (RMT) for the energy-dependent scattering 
matrix 5(e) and refer reader to review a 20 ! 21 for details 
of RMT and its relation to Green function technique. 
The diagrammatic technique based on small parameter 
1 /N <C 1 expansion are given in^ 2 . for energy- independent 
scattering matrix 5 and in 2 ^ for energy-dependent 5(e). 
We also require that the transport is only weakly non- 
linear and refer to2i24 for a discussion of highly nonlin- 
ear transport in mesoscopic samples. Here we assume 
that eV <§C TV A and treat the I — V nonlinearity only 
to order (eV) 2 , see Eq. (flj. The energy dependence 
of the 5-matrix allows us to consider non-zero tempera- 
tures T <C -E-Th, and for convenience we normalize it to 
dimensionless parameter t = 2ttT/(NA). 

The magnetic field appears via the total flux <j> through 
the dot. For simplicity the dot is assumed here to be 
(nearly) circular with radius L. As discussed in Sec. U 
the relevant flux-scale for the dot in the diffusive and the 
ballistic regime is $ c <C $o n ' 20 - Therefore for conve- 
nience we normalize magnetic flux $ through the dot to 
dimensionless flux <b as follows: 
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To describe the crossover from the diffusive to the bal- 
listic regime we use the substitute I — > 7rL/4 according 
to2£ (the numerical factor given i n 20 i 26 is not correct). 
The importance of crossover is evident if one notices that 



much of the data^ taken pertains to this intermediate 
regime. 

The geometrical capacitance C of the dot can be de- 
fined, e.g. in the constant interaction model in the Hamil- 
tonian approach 2 ^ from the solution of electrostatic prob- 
lem, if the shape of the dot is known exactly and all 
nearby conductors are taken into account. Usually such 
an exact solution is not available, but this capacitance 
can be estimated either as C ~ eL for an ungated sam- 
ple or as C ~ eL 2 /d for macroscopic gates separated 
by a distance d from the dot. (Here e is the dielectric 
constant). This capacitance defines the strength of the 
Coulomb interaction, and electronic repulsion is referred 
here as 'strong' if a typical charging energy is large com- 
pared to the level spacing, e 2 /C 3> A. 

For strong screening in the dot, r s — (fcpas) -1 = 
e 2 /(ehvF) < 1, an RPA treatment of Coulomb interac- 
tions is sufficient. Furthermore, for large dots, L clb, 
only the long distance screening is relevant. As a con- 
sequence, electrons with kinetic energy e have a well- 
defined electro-chemical potential e a = e — eV a in the 
contact a and e = e — eU(r) in the dot. For a quan- 
tum dot large compared to the Bohr radius but still 
so small that its dimensionless conductance, gdot = 
i?Th/A is much larger then conductance of the contacts, 
<?dot S> N, the potential can be taken uniform ('zero- 
mode approximation') 2 ^. The leading interactions are 
then present in the form of a Hartree electrical potential 
U which shifts the bottom of the energy band in the dot, 
see Fig. [3] and therefore modifies the 5-matrix. 

Therefore, transport depends on the Fermi- 
distributions f{e a ) and the scattering matrix 5(e). 
The current in contact a is I a = J del a {e) and for 
eV <C iVA the spectral current I a (e) can be expanded 
in powers of eV: 
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In Eq. the total current I a is expressed in terms 
of the dimensionless linear conductance at energy e, 
g<x0(s) = kr(TL a 5 a p — 1^5' t (e)l Q 5(e)) and the nonlin- 
ear conductance (measured in units of inverse energy) 
g<x0i{£) = {h/2i/ s e 3 )d 2 I a (e)/dV{3dV 11 which depends on 
the potential U. To this accuracy U needs to be known 
only up to the first order derivatives, the characteris- 
tic potentials u a = dU/dVo t 27 ^ 28 ^ 9 . The characteristic 
potentials u a are found self-consistently 30 from current 
conservation and gauge-invariance requirements. The dc 
limit ijj — > of an ac-result^i reads 



-J def'{e)trS^d E Sl a 
2iriC/v s e 2 - Jde/'(e)tr5ta e 5' 

2TtiC/v s e 2 
2iriC/v s e 2 - J de/'(e)tr5ta £ 5 = 
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We point out that this self-consistent treatment of 
Coulomb interactions has certain limitations. Remaining 
in a single-particle picture of scattering, we neglect many- 
body effects. We consider here a dot with multi-channel 
ballistic contacts, N ^> 1 when scattered wave-packets 
are delocalized. Then the effects of Fock terms are 
small < l/N compared to that of Hartree self-consistent 
potentia&2£. 

The sample-specific fluctuations of u a are dependent 
on magnetic flux (ft, and determine the field-asymmetry 
of the non-linear conductance. These derivatives are used 
to express the conductances g a /3~f (,£)'• 
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where the prime stands for energy derivative. Transport 
coefficients in general depend on magnetic field due to 
sensitivity of the matrix S to magnetic field. For ex- 
ample, field inversion transposes the scattering matrix 
S((f>) = S T (—(f>). The response to the gate voltage re- 
mains invariant under such inversion, ito(— </>) — uo((ft), 
but the response to the voltage at the contacts is asym- 
metric, u a {4>) ^ «„(- (ft). 

We point out that this 0-sensitivity of electrostatic po- 
tential U(f, (ft), which locally shifts the bottom of the en- 
ergy band, is a general feature of coherent systems. In 
the high-field limit of Quantum Hall bar a voltage applied 
to one of the contacts changes the electric potential only 
in the outgoing edge, but not in the incoming. The field 
inversion reverses the direction of edge channels, so the 
same voltage induces a potential on the opposite edge, 
U(r, $) 7^ U(f, — $). The magnetic field asymmetry in 
an edge state geometry with a Coulomb blockaded im- 
purity is the subject of Ref. HH. Even at the small fields, 
$ < $0i of interest here, the potential U(r) remains field- 
sensitive, although the field effect is not as drastic as in 
the edge state geometries mentioned above. 

Since the dot is a good conductor compared to con- 
tacts, i?Th/A ^S> N, the voltage drops mainly over the 
ballistic contacts and remains uniform inside the dot. 
This is qualitatively sketched by the bold curve in Fig. 
[3] (dashed curves correspond to the potential for chan- 
nels reflected by the ballistic constriction). Therefore 
in a quantum dot we deal with a uniform potential 
U(r,(ft) = U((ft). This is very different from an open 
diffusive sample, where r crg ~ and the voltage drops 
gradually. 

Classically the voltage through such a conducting 
nod is divided according to the widths of the contacts, 
Mela = dU c \/dV a = N a /N ~ 1, is insensitive to magnetic 
field. However, since the phases of scattering amplitudes 
between different channels are modified by the mag- 
netic field, a field-sensitive and sample-specific correction 
5u a = u a — u c \ a ~ l/N appears. Although at such small 
fields it is small compared to classical the value, its asym- 
metry with respect to magnetic field du a ((ft) ^ 8u a (—(ft) 
leads to the asymmetry in the current- voltage character- 
istic, and more specifically in the non-linear conductance 



FIG. 3: Electric potential U (x) shifts the bottom of the energy 
band. Bold curve corresponds to open channels and dashed 
(dot-dashed) correspond to (classically) reflected channels. 
Potential U varies only in the contact regions and remains 
constant inside the dot. 



9aj3i{4>) 7^ 9ap-y{~<ft)- Below, these (anti-)symmetric 
components of conductance, (G a )G s , 



^} I de(g a/3j (e,(ft)±g a p J (e,-(ft)), (7) 



are investigated in detail (for a 2-terminal dot Eq. ([7]) 
corresponds to Eq. (fTJ) for a = /3 = 7 = 1). 



III. RESULTS 

We group our results into two subsections. First in Sec. 
IIII Al we consider two-terminal dots and discuss proper- 
ties of the ensemble averaged non-linear conductances 
Qa-iQs- Since the field- asymmetry of the measured sig- 
nal is characterized by the relative strength of the com- 
ponents introduced in Eq. (0, we start by considering 
fluctuations and correlations of A = Ga/Gs, and discuss 
the dependence of A, Ga , G s on contact width and arbi- 
trary interactions. Second, in Sec. IIII Bl we generalize 
our treatment of G a to multi-terminal quantum dots and 
discuss the role of a gate voltage. 



A. Two-terminal dots 

Here we consider the (anti-) symmetric nonlinear con- 
ductance in two-terminal quantum dots in terms of the 
5-matrix and its energy derivative d £ S. To be definite, 
we set V2 — and consider derivatives with respect to 
V\ only (the gate voltage V Q is considered in Sec. IIII B|) . 
With this set of voltages, the non-linear conductances 
are G a (s) = G a (s),m- One convenient representation of 
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Qa.s is given in terms of a traceless matrix A, which is 
often used for the (dimensionless) linear conductance g 
of a quantum dot 



iV 



iV 



to separate the classical non-fluctuating part from the 
quantum contribution, the second term in Eq. ([5]). The 
projection operator 1 Q in Eq.© corresponds to the di- 
agonal matrix with a unit block for channels of the a-th 
lead and zero otherwise. Another equivalent way is to 
use the (real) injectivity d £ n a and the emissivity d e n a of 
the dot into and out of the contact a 27 ' 29 : 



d £ n a = - — : tr TL a S^d e S, 
2m 

deRa = -^—trS^ladeS. 



(9) 
(10) 



Using the matrix A or Eq. (J9j) we factorize Q a , Q s into 
products of fluctuating quantities: 



Qa,. 



tt JJdeds'f'(e)f'(e') Xl (e) X2Ms) (e') 



in) 



A 2 C/(e 2 u s ) - Jdef'(e)tT&d e S/(2m) 
Xi{e) = -^trAStAS, (12) 

X2.a(e) = ^-.tvA[d e S ) S^}=A(d e n 1 -d E fi 1 ), (13) 

X2, s (e) = ^7tr(l a St-St]Li)0 e S+i^ 
Ztxi e i/« 



A d e n 2 - d e n x 



C 



(14) 



The injectivity and emissivity are related via T, a d E n a = 
s c,9 e n Q = trS^d s S/(2m) and d e n a {<j)) = 9 E n a (-0). 
Their mesoscopic averages (d E n a ) — (d e n^) = N a /(NA) 
are denned by the coupling of the dot to the contact 
a. In contrast to their classical values, there are quan- 
tum fluctuations that render d E n a ^ d e n a and thus lead 
to non-zero r.h.s. in Eq. (fl~3|) . These averages readily 
demonstrate that (x2a) = 0, but (x2s) — CA/(v s e 2 ) + 
{N 2 -N 1 )/N. 

One can also demonstrate that the mesoscopically av- 
eraged (Q a ), (Qs) = using averaging over an ensemble 
of energy-dependent matrices (which are also symmet- 
ric in the absence of a magnetic field). Either compo- 
nent of Q is expressed by Eq. (jlip as a combination of 
fluctuating quantities, the most important being fluctu- 
ations of the numerator. One can notice from Eqs. |5J 
[T!?)) that xi(e) cx d e G(e). The latter was considered and 
proven to vanish upon ensemble averaging in Ref. l34l . 
This leads to a vanishing of the ensemble average of the 
non-linear conductance (Q a ), (Q s ) = 0. Consequently the 
magnetic field asymmetry discussed here is a signature 
of a quantum effect, similarly to e.g. the dc current gen- 
erated by an ac voltage 8 . (Interestingly, also quantum 
pumping 3 ^ exhibits a magnetic field asymmetry 3 ^ sim- 
ilar to the one discussed here.) Thus the magnitudes 



of Q a ,Qs have to be described in terms of fluctuations 
or pair-correlation functions averaged over a relevant en- 
semble of 5-matrices. 

Averaging of S(s) is performed using the diagrammatic 
technique of an expansion in l/N <C l 22 ' 23 . For ballistic 
contacts the results greatly simplify, since energy- and 
flux-dependent elements of the scattering matrix vanish, 
(Sij(e,<p)) = 0, and the non- vanishing pair correlator 
reads 

(5y(£,0)5fc ; (e',0')) = 5 ik 5jiV e - £ > + 8 a 5 jk C £ - e . , (15) 



D e — e i J N 1 + (0±0') 2 -2ni(e- e')/NA 



(16) 



The 4th-order correlator can be also expressed in terms of 
the Cooperon C e _ e < ((f), <j>') and the Diffuson T> K - F i (</>, d>') 
(higher-order correlators are presented in Refs. 1511371). 
We point out that in the Green function technique with 
random Hamiltonian the correlators (Gfj (e, (t>)Gf k (e', <fi')) 
are given by Eq. (fl"5|) (up to a normalization factor)^!. 

The denominator of Eq. (fTT|) is a self-averaging quan- 
tity, ((...) 2 ) = ((...)) 2 = (C/(C^A)) 2 , with the electro- 
chemical capacitance C M = C/(l + CA/(v s e 2 ))^ (A/v s 
stands for the mean level spacing of spin-degenerate sys- 
tem). The crossover from weak to strong interaction 
corresponds to the increase of C^/C from C^/C — > 
u s e 2 /CA fa to 1. The diagrammatic technique for 
matrices S(e) proves that the functions xi( £ i0) and 
X2a.2s(z' , 4>') are mutually uncorrelated, so that only cor- 
relators (xL'Xia) are important here. 

The simplest quantity of interest is the relative asym- 
metry, the ratio A = Q a /Qs, which strongly depends on 
the interplay between the interaction strength and the 
asymmetry of the contacts. In an experiment it may be 
desirable to maximize the magnitude of A. This quan- 
tity depends only on X2a,2s and vanishes on the average, 
(A) =0. Its correlations are non-trivial and as func- 
tions of temperature t, t' = (2tt/NA)T, T and flux 0, ft 
denned in Eq. ([J]) they are 



(A(t,ft)A(t',ft)) 



T± — T-d ± Tc 



T- 



N 2 



( C_ _ 2Ni_ 

2N X N 2 \C„ N 

00 ntx Tit' x 



2 i 



(17) 



1 

N 2 J sinh7rte sinh7rt'a; 



e x dx 



~(<t>-<t>') 



1 



A')2 ± I 



A'^2 



Note that due to interactions C^/C ^ the expression 
(|17p is not symmetric with respect to N\ «-> N% which re- 
flects the fact that we deal with a non-equilibrium trans- 
port coefficient. 

First we consider high- field limit \<f>+ft\ 3> 1 » \<j>— <f>'\, 
when the Cooperon contribution vanishes Tc *C Td so 
that T- — > T + . In case of strong interactions in a dot 
with symmetric contacts C/C M — 2N\/N — > 0, so that 
(A(t, cj))A(t' , ft)) — > 1 and Q a ,Qs could be plotted on the 
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FIG. 4: Variance of A — Ga/Gs as a function of contact asym- 
metry reaches maximum at Ni/N = C/2C M (presented for 
N — 4 and C/C^ = 1.5). As TV — > co maximum becomes 
sharper 



same scale, see Fig. [2] If either interactions become 
weak or the contacts become strongly asymmetric, the 
magnetic-field asymmetry vanishes, A — * 0. 

For strong interactions even a small asymmetry in the 
multi-channel contacts N± = N/2 ± 1 can make the sec- 
ond term in the denominator of Eq. (|17|) dominant and 
A. —> 0. The important role of contact asymmetry be- 
comes clear from Eqs. (| 131 14|) : quantum fluctuations 
6x2 ~ min {1, y/l/t} /N of numerator and denominator 
of A ~ Xialxis occur on the background of (x2a) = 
and \X2s) = C/C/i ~ 2N\/N. The numerator grows pro- 
portionally to 8x2- In contrast for the denominator the 
fluctuations can be totally neglected if its classical av- 
erage value is large (x2s) ^S> 8x2 and therefore A -C 1. 

In an experiment the number of channels can be ad- 
justed to maximize asymmetry A. In a weakly interact- 
ing dot with C/Cfj, > 2 the maximum of A 2 is reached 
at N x /N = C/2(C - C p ), but it is still small A < 1 due 
to the weak interactions. If the interaction is relatively 
strong, C/C M < 2 then at Ni/N — C jlC^ the maximum 
A = 1 is achieved, see Fig. @] It is this contact width 
that minimizes the magnitude of (x2s)- 

In a realistic quantum dot with symmetric contacts^ 
the interaction is reasonably strong and at intermedi- 
ate magnetic fields \(j> ± <f)'\ < 1 we have A 2 = T-j T+. 
However, it is important to emphasize that if a dot be- 
comes smaller, the ratio e 2 /CA (oc L/as for a dot with- 
out gates) can diminish as well. Therefore one could 
expect that interactions are not necessarily strong and 
at 1 < C/Cp < 2 it would be important to make the 
contacts asymmetric to clearly observe a magnetic-field 
asymmetry. The peak in the Fig. [4] becomes narrower 
as N — * 00 and the choice of contact width such that 
Ni/N sa C/2C^ becomes crucial. This behavior is a con- 
sequence of the scaling of the numerator of Eq. (|17p : the 



width of this peak as a function of N\/N is proportional 
to fluctuations 8x2 discussed above. Although experi- 
ments are often performed with symmetric contacts to 
avoid parasitic effects of an asymmetric circuit, our re- 
sults suggest that using contact asymmetry can enhance 
the magnetic-field asymmetry. 

Fluctuations of A can not be found in closed form, 
since T± are complicated functions of t and cj>, see Eq. 
(fTT|) . At low temperatures t -C 1 the temperature- 
dependence of T± can be neglected. In the regime t > 1, 
more relevant for experiment^ (see also discussion in Sec. 
lIVp . the fluctuations of A for a dot with symmetric con- 
tacts Ni = N2 can be expressed in terms of linear con- 
ductance g(4>) and its universal fluctuations var g: 



var A 



var A 



var g 



(18) 



(2 + ^ 2 )var 5 + 2(l + </» 2 )(^--l) g* 

1 0>1. (19) 



l + 2( 5 7var gXC/C^-l) 2 



The high-field limit <fi 3> 1 of var A is universal and 
Eq. (fT9")) is valid for arbitrary temperatures t, as long as 
T <C -&rh- However, this value is temperature-dependent 
due to temperature dependence of UCF. As Eqs. (|18I19|) 
demonstrate, measurements of A could serve as a tool to 
find interaction strength Ca/C. Recently an experiment 
in micrometer-sized Aharonov-Bohm (AB) rings^ found 
A ~ 0.3 analyzing the amplitude of the AB-oscillations in 
the non-linear conductance. This asymmetry was found 
to be rather insensitive to variation of sample's conduc- 
tance j > 1, which means that the 2nd term in the 
denominator of Eq. (|19p is not large and interactions are 
rather strong, C^/C w 0.9. 

Since the asymmetry parameter A is not sensitive to 
the fluctuations of the density of states, it probes in- 
teractions more directly than Q a . The symmetrized Q s 
is measured ir>i^ on the background of UCF in the lin- 
ear conductance, so that the interaction strength could 
be evaluated only from the amplitude of Q a . Similar to^ 
measurement of A in quantum dots would provide a more 
precise measure of interactions then those of Q a , Q s , which 
we consider below. 

For pair correlators of Q a and Q s all fluctuating quan- 
tities of Eq. (TTTj) should be taken into account: 



(Ga(t,4>)Ga(t', <!>')) 

(g s (t,,t>)g s (t>,4>>)) 



N 2 f _C Ml 

VC M N 



A C 



2N!N 2 
2 



~lX7G~ 



T. 



+ • 



(20) 



Some properties of these correlators at t = t', <f> = <f>' were 
discussed in Ref. [ll|. Equation (|2H)) allows one to easily 
evaluate fluctuations in the non-linear conductance as a 
function of temperature and flux. For example, at suf- 
ficiently high flux difference, correlation functions decay 
strongly. Similarly to 'magnetic fingerprints' in the linear 
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conductance^, the nonlinear conductance is randomized 
beyond \4> — <j>'\ ~ 1. 

We note some similarities between Eqs. <\17\ and 
Eq. (|2"0")) . Despite these similarities, an important differ- 
ence is that correlations of A = G a /Qs (see (HZ])) are not 
given by a simple ratio of its averaged numerator and 
denominator of Eq. {2D]), (A 2 ) ^ {Ql)/{Q 2 ). In contrast, 
to leading order in 1/N <C 1 Q a and Q s are uncorrelated. 

Weakly interacting limit. In the limit of weak inter- 
actions C^/C — > and the magnetic field asymmetry 
in Eq. (|20[) vanishes. However, even without interac- 
tions non-linear conductance exists^. The energy scale 
eV c ~ ft/id was argued to be the only energy scale for 
the crossover from linear to non-linear transport. In the 
open diffusive samples considered in&£ eV c ~ i?Th since 

7"d ^ 7"erg- 

Equation (|20p helps us to relate the correlations of non- 
linear conductances at arbitrary t, <p and t', <f>' in weakly 
interacting dots to those of the linear conductance. In 
particular, we find that in this limit the ratio of corre- 
lation functions of linear and non-linear conductances is 
constant at fixed Td. Particularly, for UCF and fluctua- 
tions of G 7 we find a universal relation 



var g = 2{h/T^) 2 var Q s 



(21) 



insensitive to temperature t and magnetic flux <f>. This 
relation illustrates the qualitative arguments of&£ that 
transport coefficients in the non-interacting system scale 
with h/r d = NA/2n. 

In the crossover from weak to strong interactions a 
typical amplitude of the fluctuations in Q s changes from 
SQs ~ l/(NA),Sg a « 0, see Eq. to a value 

Sg s ,8g a ~ 1/(AT 2 A) . min {1, y/TJt}. Therefore, in the 
experiment, the interaction induced scaling of a non- 
linear response SQ a ~ 5Q S proportional to 1/N 2 can be 
observed. In contrast for a non-interacting system we 
predict (see Eq. (|2"TjO that Q a = and Q s is proportional 
to 1/N. Indeed, in the experiment by Angers et al*£, 
this observation allows one to evaluate the interaction 
strength in Aharonov-Bohm rings, C^/C w 0.9. 



B. Multi-terminal dot 

We point out that similarly to two-terminal dots one 
can consider multi-terminal dots as well. Since Q a is very 
sensitive to interactions and is often easier to find ex- 
perimentally, we concentrate on a similar anti-symmetric 
component Q a ,aaa and on a mixed non-linear conduc- 
tance G a ,af3Q, a derivative with respect to the gate voltage 
Vo and contact voltage Vp. 

One can demonstrate that in two-terminal dots the 
non-linearity with respect to Vq is absent. Indeed, Q a /3Q in 
Eq. ([6]) is only due to Uo(4>)g' a3 (<p). The derivative ug is 
always field-symmetric, see Eq. In two-terminal dots 
9' a p{4>) is symmetric too, g' a3 {<j>) = g' a p{-4>)- A s a con- 
sequence Gaps is even in magnetic field for two-terminal 



dot (the classical effect of Vo in two-terminal samples was 
considered ir*2£). However, such asymmetry appears in a 
multi-terminal set-up, since g' a g{4>) ^ g'afii"^)- 
We find for the variances of Q a ,aaa and G a , a /30i 



Var Ga.aaa — 
Var Ga^a/30 = 



fcrCpV N*(N - N a f 
A C 



2 N 2 N 2 
A 4 



A 6 
T-. 



T+T-,(22) 
(23) 



Equation (|23p is of particular importance since it can 
serve to find C M in an experiment. This result does not 
contain an apriori unknown C. It is useful to notice that 
for strongly interacting electrons the effect of gate volt- 
ages is weak (C — > 0) and therefore voltage difference 
between gates and the dot can be rather large. There- 
fore, the non-linear conductance becomes insensitive to 
the gate voltage and scales with C M w C — > 0. However, 
for weakly interacting electrons C M — > v s e 2 /A, the inter- 
nal potential closely follows the gate voltage and so the 
effect of gates becomes measurable. 



IV. DISCUSSION 

The results we have obtained can now be compared 
with results for the linear conductance. From such a com- 
parison we can assert that it is indeed the weakly non- 
linear transport regime which should be used to extract 
information on the interaction strength rather than the 
linear conductance. At finite temperature t — 2irT/NA, 
in open quantum dots with many ballistic channels, N S> 
1, the amplitude of the random linear conductance fluctu- 
ations calculated in the strong interaction limit e 2 /CA ^> 
1 depends on interactions as Sg ~ — t 2 , for i « 1 and 
Sg — -1/Ai 3/2 for t > l 5 . On the other hand, random 
fluctuations of the non-linear conductances Q a , Q s in this 
limit e 2 /CA » 1 are oc 1/(A 2 A) min {1,1/t}- For a 
voltage of the order of eV oc NA, the fluctuating con- 
tribution of the non-linear component to the current is 
oc Q(eV/A) 2 and can be compared with the fluctuating 
linear contribution Sg eV/A. We conclude that for volt- 
ages in the interval min {Nt 2 ,l/Vt} < eV/NA < 1 
the transport is still weakly non-linear, but the quantum 
contribution due to Coulomb interaction to non-linear 
transport is stronger then the linear contribution. This 
estimate demonstrates that in dc-transport the effect of 
Coulomb interactions in open mesoscopic conductors is 
best investigated in the non-linear regime. 

Another important issue is the sensitivity to small 
magnetic flux, <f> <C 1, when I(cf>) in a two-terminal dot 
becomes an asymmetric function due to Coulomb inter- 
action. In the experiment t > 0.5^, so temperature 
should be taken into account. Although we do not con- 
sider inelastic scattering due to T ^ 0, the non-linear 
conductance G a {T) is temperature dependent since scat- 
tering is energy dependent. Of particular interest for the 
experiment^ is the quantity gs- = (2is s e 3 /h)VQ a in the 
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limit where Ga{4>) is still linear in flux, </> < 1, see Fig. [2] 
We use Eq. ^ and find 



rms G a {4>) 



A C \ N 10 



1/2 



4</> 



TrH 4 x 2 y 2 (l + y)e- x - y dxdy\ 
2 sinh 2 nty sinh 2 ntx 



1/2 

) (24) 



The data obtained at T = AfieV « A = 7/ieV were fit- 
ted to a combination of the zero-temperature results*^ 
corresponding to a rms Q a given by the first line of (|2"4"|) 
with 40 substituted by $/$oH, 

It is noteworthy that a treatment of the crossover 
leads instead to Eq. (f24|) with a non-linear conductance 
rms Q a oc 0, where 4> oc $/$c is defined by Eq. ([2|), and 
not ^/^o- Indeed, the flux dependence of G a {<&) starts to 
saturate at a scale which is parametrically smaller then 
$o- This is a consequence of the existence of two rele- 
vant time-scales, and r crg . Thus an attempt to fit the 
data to a theory of open diffusive samples 10 by substi- 
tuting r org — > Td fails to catch the difference between $o 
and <& c oc $o( T erg/ T d) • Being an order-of-magnitude 
estimate to ballistic quantum dot results, the substitu- 
tion <& c — > $o misses parametric difference between these 
fields due to the zero-dimensional physics of a chaotic 
dot. This is especially important if the dependence of 
the effect on N is considered, and rms Q a oc N~ 5 / 2 at 



T ^ 0. 

We also point out that the temperature effect should 
have appeared even at such low temperatures as T — 
4/xeV. Indeed, the last factor of Eq. (|24|) is equal to 1 
for t <C 1 and \f2-K jYlt at t 3> 1. Because of exponential 
suppression of the integrands at x, y ~ 1/ (7r£), for t > 0.5 
the high temperature asymptotic is a good approxima- 
tion. Since <I> C oc V~N, the nonlinear conductance as a 
function of N and T behaves as rms £/ Q oc N~ 3 ^ 2 T~ 1 . 

To summarize, we discussed the symmetric and 
anti-symmetric components of non-linear conductance 
through chaotic two-terminal and multi-terminal quan- 
tum dots. Their correlations were found to strongly de- 
pend on the interplay between the asymmetry of ballistic 
contacts and the strength of Coulomb interactions. We 
discussed the effect of a variation of the gate-voltage on 
the anti-symmetric component of weakly non-linear con- 
ductance through multi-terminal dots. We investigated 
the conditions on sample contacts, magnetic field and 
temperature under which the (relative) magnetic field 
asymmetry is maximal and thus most easy to detect. 
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ported by the Swiss National Science Foundation, the 
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